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Abstract 


Formulae  arc  recommended  for  calculating  the  extent  of  the  spread 
of  oil  slicks  on  water  as  a  function  of  time.  They  are  based  on 
empirical  measurements  of  spreading  rates  and  analytical  and  theoretical 
studies  of  the  physical  processes  which  accelerate  or  retard  the  spread 
of  a  film.  Both  one-dimensional  and  two-dimensional  (axisymmetric) 
slicks  are  treated.  Comparisons  of  the  recommended  formulae  are  made 
with  the  limited  number  of  field  observations,  both  for  the  rate  of 
spread  and  the  maximum  slick  size. 
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I.  Introduction 


This  paper  reviews  our  current  understanding  of  the  physical 
processes  which  initially  cause  and  ultimately  terminate  the  spread  of 
oil  (or  other  immiscible  fluids)  on  the  surface  of  water.  We  are 
principally  concerned  with  the  spread  of  large  volumes  of  oil,  such 
as  might  be  encountered  in  spills  from  ships  or  oil  wells,  and  which 
cannot  be  reproduced  easily  in  a  laboratory  at  full  scale.  Our 
approach  is  to  consider  some  simple  cases  of  oil  spread,  which  will 
not  likely  be  duplicated  precisely  in  practice,  but  for  which  a 
theoretical  or  semi-empirical  description  can  be  found,  especially 
through  use  of  properly  designed  laboratory  experiments  which  simulate 
full  scale  spreading  phenomena.  Based  on  this  understanding,  a 
correlation  of  field  observations  is  made  as  a  test  of  the  suitability 
and  accuracy  of  these  predictions,  and  empirical  formulae  for 
estimating  spreading  rates  are  recommended. 

The  physically  most  important  assumption  underlying  our  analysis, 
which  is  most  likely  to  be  violated  in  any  real  incident  of  a  spill, 
is  the  absence  of  any  effects  of  wind,  tidal  currents  and  waves.  We 
would  expect  that  the  drifting  motion  caused  by  winds  and  tidal 
currents  would  simply  be  superimposed  on  the  spreading  motion  to  be 
experienced  on  calm,  stationary  water,  since  these  latter  motions  are 
confined  to  a  layer  near  the  surface  which  is  relatively  thin  com¬ 
pared  with  that  subject  to  wind  friction  effects  or  tidal,  motion. 

It  is  more  likely  that  wind  and  tidal  current  will  produce  relative 
shearing  motion  in  the  plane  of  the  water  surface,  deforming  the 
shape  of  the  spreading  slick  from  those  simple  shapes  expected  in 
calm  water.  Such  distortion  is  commonly  observed,  and  is  most  likely 
to  limit  the  usefulness  of  the  spreading  laws  which  we  propose. 

These  effects  are  very  difficult  to  predict  or  even  describe,  and 
there  is  little  empirical  evidence  on  which  to  base  an  estimate  of 
their  importance. 

The  first  order  effects  of  surface  waves,  on  the  other  hand, 
can  be  shown  to  be  negligible.  Because  of  their  periodic  nature, 
waves  produce  oscillatory  forces  having  zero  mean  value  and  which 
therefore  do  not  affect  the  spreading  motion  which  proceeds  on  a 
much  longer  time  scale  than  the  usual  wave  periods.  Of  course,  there 
are  non-linear  effects  of  waves  which  are  probably  not  separately 
distinguishable  from  those  associated  with  winds  and  tidal  currents, 
and  are  equally  difficult  to  predict. 

The  use  of  laboratory  scale  experiments  to  establish  empirical 
spreading,  laws  and  to  check  theoretical  predictions  has  been  discussed 
elsewhere1.  We  shall  make  use  of  such  experimental  evidence  to  pro¬ 
vide  the  best  estimate  of  spreading  rates,  accepting  the  asserted 
validity  of  the  scaling  of  these  experiments  to  full  size.  The  basis 
of  such  scaling  is  a  well  understood  aspect  of  fluid  mechanics,  and 
will  not  be  further  discussed  here. 


Accurate  field  observations  of  the  spread  of  oil  slicks  are 
very  rare.  We  have  tried  to  include  all  measurements  which  have  been 
published,  but  in  most  cases  we  have  been  forced  to  assume  additional 
information,  such  as  spreading  coefficients,  when  comparing  these 
observations  with  a  theory.  Given  the  inaccuracy  of  the  observations, 
these  assumptions  are  of  no  great  significance,  but  only  serve  to 
emphasize  the  scarcity  and  crudity  of  the  observations.  A  major  goal 
of  our  proposed  correlations  is  to  permit  the  comparison  with  future 
(and  hopefully  more  accurate)  observations. 

II.  Spreading  and  Retarding  Forces 

Although  the  force  of  gravity  acts  downward,  it  causes  a  side- 
wise  spreading  motion  of  a  floating  oil  film  by  creating  an 
unbalanced  pressure  distribution  in  the  pool  of  oil  and  the  surround¬ 
ing  water.  This  force  on  an  element  of  oil  film  acts  in  the 
direction  of  decreasing  film  thickness  and  is  proportional  to  the 
thickness,  its  gradient,  and  the  difference  in  density  between  oil 
and  water.  (See  Fig.  1.)  As  the  oil  film  spreads  and  becomes 
thinner,  the  gravity  force  diminishes. 

At  the  front  edge  of  the  expanding  slick  an  unbalance  exists 
between  the  surface  tension  at  the  water-air  interface  and  the  sum 
of  surface  tensions  at  the  oil-air  and  oil-water  interfaces.  The 
net  difference,  called  the  spreading  coefficient,  is  a  force  which 
acts  at  the  edge  of  the  film,  pulling  it  outwards.  This  spreading 
force  does  not  depend  upon  the  film  thickness  as  does  the  gravity 
force,  and  will  not  decrease  as  the  oil  film  thins  out  (unless  the 
chemical  properties  change  through  aging) .  Eventually  the  surface 
tension  force  will  predominate  as  the  spreading  force. 

These  spreading  forces  are  counterbalanced  by  the  inertia  of 
the  oil  film  and  of  the  thin  boundary  layer  of  water  below  it  which 
is  dragged  along  by  friction  (see  Fig.  1).  The  inertia  of  an 
element  of  the  oil  layer  decreased  with  its  thickness  as  time 
progresses  and  the  film  spreads,  but  the  inertia  of  the  viscous 
layer  of  water  below  the  oil  increases  with  time  as  its  thickness 
grows.  Consequently,  the  viscous  retardation  will  eventually  out¬ 
weigh  the  inertial  resistance  cf  the  oil  layer  itself. 

It  is  also  informative  to  consider  these  effects  from  the  point 
of  view  of  an  energy  balance.  A  pool  of  oil  floating  on  water 
possesses  a  greater  potential  energy  thar  the  water  it  displaces, 
in  proportion  to  its  thickness.  As  it  spreads  and  its  thickness 
decreases,  there  is  a  loss  of  potential  energy.  Also,  as  air/vater 
surface  is  replaced  by  an  oil  film,  the  surface  energy  per  unit 
area  (which  has  the  same  physical  value  as  the  interfacial  tension) 
is  reduced  by  an  amount  equal  to  the  spreading  coefficient.  Thus 
both  surface  energy  and  potential  energy  are  decreased  as  the  slick 
spreads.  This  energy  is  converted  either  into  heat  by  viscous 
dissipation  in  the  water  beneath  the  slick  or  into  the  energy  of 
gravity  surface  waves  which  propagate  awav  from  the  expanding  oil 
pool.  In  other  words,  each  spreading  force  is  associated  with  an 
energy-producing  process  and  each  regarding  force  with  an  energy- 
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dissipating  process. 


It  is  thus  clear  that  the  spread  of  an  oil  film  will  pass  through 
several  stages  as  time  progresses,  in  each  of2whlch  one  spreading 
force  will  be  balanced  by  one  retarding  force  .  Although  there  are 
four  such  possible  combinations,  for  large  scale  slicks  only  three 
regimes  are  important:  (i)  the  gravity-inertia  regime  (called 
"inertial  spread"),  (li)  the  gravity-viscous  regime  (called  "viscous 
spread")  and  the  surface  tension-viscous  regime  (called  "surface 
tension  spread").  As  time  progresses,  a  large  spill  will  pass  through 
these  three  regimes  in  succession.  A  very  small  spill  (a  few  liters, 
say)  will  almost  form  the  start  behave  as  a  surface  tension  spread. 

The  spreading  laws  for  each  regime  have  been  determined,  to 
within  an  unknown  constant,  for  each  regime  and  for  the  ^a^es  of  a 
one-dimensional  and  two-dimensional  (axisymmetric)  slick  ’  .  These 
laws  give  the  linear  extent  of  the  slick  (length  Z  or  radius  r)  as 
a  function  of  the  time  t  since  the  oil  was  released  at  the  origin 
of  the  spread,  the  volume  of  the  oil  spill  and  the  physical  properties 
of  the  oil  and  water.  These  spreading  laws  are  given  in  Table  I, 
and  the  undetermined  proportionality  coefficients  are  denoted  by  the 
symbol  k. 

III.  evaluation  of  Spreading  Laws 

The  proportionality  constants  1<  can  he  determined  from  laboratory 
experiments  or  from  a  suitable  detailed  hvdrodvnamic  theorv  of  the 
spreading  motion  in  each  regime,  ^o^far,  only  one-dimensional  spread¬ 
ing  experiments  havn^bgen  reported  *  ,  and  there  have  been  advanced 
conflicting  theories  ’  for  the  inertial  spreading  regime.  We  suggest 
below  (and  in  Table  II)  best  values  for  these  coefficients,  based 
upon  published  experimental  data  and  our  own  (unpublished)  theoretical 
analysis  and  extrapolation  of  the  empirical  data.  We  discuss  below 
each  entry  separately. 

One-dimensional  inertial  spread.  Here  we  us£  the  experimental 
value  of  k.  ■  1.5  determined  by  Hoult  and  Suchon  (see  Fig.  2). 

Their  theoretical  value  (k.  ■  3)  is  clearly  in  disagreement  with 

the  experiments.  An  alternative  theoretical  solution  has  been 
given, hv  Fannelop  and  Waldman  ,  from  which  k.  is  found  to  be 
3/10/.J  ■  1.39.  We  believe  that  the  correcctheoretical  value  is 

3/ 7 A  ■  1.57,  for  reasons  which  we  shall  not  elaborate  on  here. 

This  latter  theoretical  estimate  is  certainly  very  close  to  the 
observed  ’*alues  in  the  laboratory  experiments. 

One-dimensional  viscous  sgread.  The  empirical  value  of  k^v  “  1.5 
determined  by  Hoult  and  Suchon'*' (see  Fig.  2)  is  recommended. 

One-dlmensiongl  surface  tension  spread.  The  experiments  of 
Garrett  and  Barger"'(see  Fig.  3)  and  Lee"'  (see  Fig.  A)  both  support 
a  value  of  *  1.33. 
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Two-dimensional  (axlsymmetric)  inertial  spread.  Since  there 
are  no  experimental  results  available,  we  recommend  the  theoretical 
value  of  k„  ■  2/(3n)1  ■  1.14  as  determined  hy  the  same  analvsis 

leading  to  the  one-dimensional  value  quoted  above,  and  which  agreed 
closely  with  the  corresponding  experiments.  This  value  is  also  given 
bv  Fannelop  and  Waldman  . 

Two-dimensional  viscous  spread.  Again,  no  experiments  have 
been  reported.  A  boundary  layer  theory  developed  by  Moult  and 
Suchon  possesses  no  unique  solutions  and  hence  yields  no  definite 
values  for  k.  or  *s  our  belief  that  the  proper  solution 

can  only  be  determined  theoretically  hv  solving  the  complicated 
flow  at  the  leading  edge  of  the  slick.  However,  we  suggest  that  an 
estimate  of  k2  can  he  made  in  the  following  manner.  If  we  select 
the  one-dimensional  theoretical  solution  which  leads  to  the  observed 
value  of  and  then  hypothesize  that  the  two-dimensional  solution 

should  have  the  same  boundary  values,  we  can  then  determine  the 
value  of  k^v  from  this  latter  solution.  Our  justification  for  such 
a  procedure  is  the  supposition  that  the  flow  near  the  leading  edge 
of  the  slick  is  the  same  for  the  two-dimensional  as  for  the  one- 
dimensional  case,  and  that  the  dimensionless  boundary  values  of  the 
theoretical  solutions,  which  are  determined  by  this  flow,  should 
also  be  identical.  Using  this  procedure,  we  have  found  the  value 
shown  in  Table  II. 

Two-dimensional  surface  tension  spread.  We  have  used  the  same 
procedure  as  that  outlined  in  the  preceding  paragraph  to  estimate 
the  value  of  k^t,  shown  in  lj)ble  II,  since  there  are  no  experiments 
available.  According  to  Fay  ,  the  maximum  observed  spread  in  field 

observations  would  correspond  to  k£t  *  lO/ir"  ■=  5.7.  While  this  is 
larger  than  the  corresponding  value  recommended  in  Table  II,  it  is 
most  likely  uncertain  by  a  factor  of  two  because  of  the  difficulties 
of  making  observations  and  the  imperfections  of  the  field  experi¬ 
ments.  A  comparison  of  the  theoretical  spread  area  with  observed 
values  is  shown  in  Fig.  5. 

IV.  The  Termination  of  Spreading 

It  h^s^already  been  noted  that,  after  some  time,  slicks  cease 
to  spread^’  .  In  almost  all  cases,  the  final  f^lm  thickness  is  _ ^ 
much  greater  that)  that  of  a  monomolecular  layer  ,  being  about  10 
to  10  cm.  Fav  has  suggested  that  the  cessation  of  spread  is 
caused  bv  the  evaporation  of  some  oil  fractions  which  reduces  the 
spreading  coefficient  to  zero.  His  estimate  of  slick  size  for 
which  this  evaporation  (limited  hy  diffusion  through  the  oil  layer) 
would  be  appreciable,  was  an  order  of  magnitude  smaller  than  the 
observed  values. 

We  propose  here  a  modified  version  of  this  theory.  We  believe 
that  the  spreading  coefficient  is  reduced  by  an  increase  in  the 
water-oil  interfacial  tension  brought  about  by  the  dissolving  of 
oil  fractions  in  the  water  layer  underneath  the  oil  film.  The  volume 
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of  oil  which  can  be  dissolved  In  this  layer  (per  unit  area  of 
oil/water  Interface)  would  then  be  s(Dt)1/  ,  where  s  Is  the 
solubility  of  the  significant  oil  fractions  In  water.  As  a  con¬ 
sequence,  th^.grevious  estimate  of  Fay  would  be  increased  by  a 
factor  of  s  ,  or  a  factor  of  about  ten  for  s  ■  10  ,  a  reasonable 

value.  As  a  consequence,  the  maximum  area  A  of  the  slick  would  be¬ 
come. 


A 


k  /q2v6  *1/8 

a'  2  _3  f>} 
p  VD  s 


(1) 


In  which  ka  Is  an  undetermined  constant  of  order  unity. 

Because  of  the  uncertainties  In  s  and  a  in  the  field  observa¬ 
tions,  and  the  lack  of  laboratory  data,  It  Is  proposed  that  the 
maximum  slick  area  be  related  to  the  volume  of  the  spill  by  the 
dimensional  formula. 


A(m2)  -  105[V(m3)]3/A 


(2) 


This  Is  compared  In  Fig.  6^with  field  observations  recently  sum¬ 
marized  by  Allen  and  Estes  .  Equation  (l)_woul«;  have  the  val^e 
given  by  Eq.  (2)  if  a  -  10  dyne/cm,  D  ■  10  5  cm  /sec,  s  -  10  , 

and  k  ■  1. 
a 


*111  J»ll|  4  Hifii  I  IWliffWmB’ 
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Table  I 


Spreading  Laws  for  Oil  Slicks 


One-dimensional  Axisymmetric 


Inertial 

t  -  ku(4gAt2)1/3 

r  -  k2i(AgVt2)1/4 

Viscous 

l  -  k^UgkV'V'2)1/* 

r  -  k2v(AgV2t3/2/v1/2 

Surface  tension 

l  -  klt(o2t3/p2v)1M 

r  -  k2t(a2t3/o2v)1/4 

Table  II 

Spreading  Law  Coefficients 

One-dimensional 

Axisymmetric 

Inertial 

m 

• 

H 

a 

•H 

•Ad 

1.14 

Viscous 

kl,  ■  l-3 

1.45 

Surface  tension 

k,t  -  1.33 

2.30 
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List  of  Symbols 


A  Volume  of  oil  per  unti  length  normal  to  x 
g  Acceleration  of  gravity 

h  Thickness  of  oil  film 

k  Proportionality  constant 
£  Length  of  one-dimensional  oil  slick 
r  Maximum  radius  of  axlsynsnetrlc  oil  slick 
t  Time  since  initiation  of  spread 
u  Spreading  velocity  of  oil  film 
V  Volume  of  oil  in  axisymmetric  spread 
x  Dimension  in  direction  of  one-dimensional  spread 
6  Thickness  of  viscous  boundary  layer  in  the  water  underneath 
the  oil  film 

o  Spreading  coefficient  or  interfacial  tension  (with  subscript) 
v  Kinematic  viscosity  of  water 
V>  Absolute  viscosity  of  water 

p  Density  of  water 

A  Ratio  of  density  difference  between  water  and  oil  to  density 
of  water 

Subscripts 

1  One-dimensional  spread 

2  Two-dimensional  (uxlsymmetric)  spread 
a  Maximum  area 

1  Inertialspread 
t  Surface  tension  spread 
v  Viscous  spread 
ow  Oil /water 
aw  Air/water 
oa  Oil/air 
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Physical 


List  of  Captions 


Hie  four  forces  which  act  on  an  oil  film  (see  list  of 
symbols) . 

Experiments  showing  the  transition  from  inertial  to 
viscous  spread  for  a  one-dimensional  flow  . 

Measurements  of  spreading  velocity  versus  slick  length 
for  ovie-dimensional  surface  tension  spreading  experi¬ 
ments'*.  - 

Lee's  experiments'3  on  one-dimensional  surface  tension 

spreading.  Solid  line  corresponds  to  k.  *  1.33. 

A  comparison  of  the  theoretical  axisymmetric  slick  ar^a^ 
(for  surface  tension  spread)  with  observed  values  *  . 
Solid  line  corresponds  to  the  value  of  k_  shown  in 
Table  II  and  a  spreading  coefficient  of  30  dyne/cm. 

Maximum  slick  area  as  a  function  of  volume.  Eq.  (2) 
compared  with  observations  taken  from  Ref.  7. 


Processes  in  the  Spread  of  Oil  on  a  Water  Surface  by  James  A.  Fay 
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Experiments  showing  the  transition  from  inertial 
to  viscous  spread  for  a  one-dimensional  flow 
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Fig.  5  A  comparison  of  the  theoretical  axisymmetric  slick  area 
(for  surface  tension  spread)  with  observed  values 
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Fig.  6  Maximum  slick  area  as  a  function  of  volume 
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